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ABSTRACT 

Derivation  and  methods  for  the  proof  of  existence  of 
singular  lines  on  the  parameter  plane  is  presented.   A  method 
is  derived  for  producing  singular  line  systems.   Several  sys- 
tems are  designed  and  evaluated  using  this  technique. 
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I .   INTRODUCTION 

Methods  for  studying  the  effect  of  variations  in  param- 
eters of  a  multivariable ,  multiloop  control  system  have  been 
developed  by  Mitrovic,  Siljak,  Thaler,  et  al  (1),  (2),  (3). 

Mitrovic's  method  consists  of  specifying  as  variables 
the  two  lowest  order  coefficients  of  the  characteristic  equa- 
tion.  Using  the  characteristic  equation  as  a  mapping  func- 
tion, constant  zeta,  omega,  and  sigma  curves  are  transformed 
from  the  s-plane  into  the  variable  coefficient  plane.   This 
technique  permits  the  adjustment  of  these  variables  so  that 
the  roots  of  the  characteristic  equation  may  be  set  at  de- 
sired locations. 

In  1964,  Siljac  extended  Mitrovic's  Method  into  the 
Coefficient  Plane  Method  (2) .   This  technique  permits  the 
utilization  of  arbitrary  pairs  of  coefficients  as  variables. 

The  Parameter  Plane  Method  is  an  extension  of  the  above 
methods,  which  permits  the  adjustment  of  system  parameters 
which  are  located  in  more  than  two  coefficients  of  the  char- 
acteristic equation  (5) ,  and  (6) .   The  Parameter  Plane  Method 
permits  the  analysis  of  complex  systems  which  mav  be  either 
linear  or  nonlinear.   However,  the  complexity  of  this  tech- 
nique requires  the  computation  of  the  curves  by  a  computer. 

In  1965,  Nutting  developed  a  computer  program  (PARAM-A) 
which  calculates  and  plots  the  Parameter  Plane  curves  for  con- 
stant zeta,  natural  frequency  and  sigma  (7).   In  general, 
analysis  of  the  resulting  curves  results  in  the  location  of 
all  of  the  roots  of  a  characteristic  equation  for  a  given 


parameter  (a, 3)  pair.   However  in  1966,  Bowie  discovered  sys- 
tems for  which  the  existing  Parameter  Plane  technicrues  did 
not  locate  all  of  the  roots  (8).   Bowie  designated  these  sys- 
tems as  singular  line  systems. 

Section  III  describes  the  mathematical  basis  of  sinqular 
lines  and  methods  for  proving  their  existence.   A  method  is 
also  developed  for  obtaining  singular  lines.   Section  IV  con- 
tains examples  of  singular  line  systems.   The  application  and 
analysis  of  singular  line  techniques  to  specific  problem 
areas  is  contained  in  Section  V.   Section  VI  concludes  this 
report  with  comments  on  existing  singular  line  systems  and 
recommendations  for  further  investigation  of  singular  line 
systems . 


II.   THE  PARAMETER  PLANE 

In  this  section  the  derivation  of  the  parameter  plane 
equations  for  the  linear  case  will  be  derived  (9)  and  gen- 
eral rules  for  mapping  points  from  the  s-plane  into  the 
parameter  plane  will  be  discussed  (3). 

A.   PARAMETER  PLANE  EQUATIONS 

Consider  the  characteristic  equation: 

n     , 
F(s)  =  I      a,  sK  =  0  (2-1) 

k=0   K 


where : 


the  coefficients  a,  (k=0 ,1 , . . . ,n)  are  real 
functions  of  two  system  parameters  desig- 
nated as  a  and  3  and  s  is  the  complex  var- 
iable.  Let  s  be  written  as: 


s  =  -£oj   +  jo)  /1-C2  (2-2) 

n    J  n 

where  oo   is  the  undamped  natural  frequency  and  £  is  the 

damping  ratio.   Letting 

6  =  cos_1(-C)  (2-3) 

and  substituting  into  the  equation  for  s  gives: 

s  =  w  (cose  +  j  sine)  =  u  e3       .  (2-4) 

n  n 

Then  s   may  be  written  as : 

s   =  oo  e-1    =oj  (cos  ke  +  jsin  k6)  .      (2-5) 
n        n 


Defining : 

T.(-c)     =    cos    k9    =    cos [k    cos"    (-?) ] 

tt    t    r\    -    sin  k9   _   sin  [k   cos      (-Q  ]  ,„    c. 

uk("?)    =  -iIHe iTHe •  (2"6) 

And   substituting    into    the    equation    for   s    results    in: 

sk   =   %[Tk(-C)    +    j/l-S2    Uk(-C)]     .  (2-7) 

It  may  be  shown  that: 

Tk(-0  =  (-l)kTk(c) 

Uk(-C)  =  (-l)kUk(?)  (2-8) 

where  T,  (c)  and  uV.(C)  are  Chebyshev  functions  of  the  first 
and  second  kinds  respectively  and  are  given  by  the  recursion 
formulae : 

Tk+i(?)-2Tk(c»  +  Tk-i(c)  =  ° 


where 


Uk+l(C)"2Uk(C)  +  Uk-l(c)  =  °  (2"9) 


T0(c)  =  1     UQ(c)  =  0 


T]_(c)  =  C    ux(c)  =  -1   • 

Substituting  equation  (2-9)  into  equation  (2-1)  and  equa- 
ting the  real  and  imaginary  parts  independently  to  zero  re- 
sults in: 

I      avo)k(-l)ETv(?)  =  0 


k=0   k  n 


I      a  k<(-D*+XU  .(C)  =  0   . 
k=0  K  n       * 


(2-10) 
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From  equation  (2-9)  the  interrelation  between  the  Chebyshev 
functions  is  obtained: 

Tk(c)  =  CUk  -  Uk_x(c)  .  (2-11) 

Substituting  into  equation  (2-10) : 

j  (-D VX-i^  =  o 

k=0 

(2-12) 

I       (-1)\A  (U  =  0  . 
k=0       K  n  K 

Assume  for  the  linear  case  that  the  coefficients  a,  are  lin- 

k 

ear  combinations  of  the  system  parameters  a  and  3  and  are  of 
the  form: 

ak  =  bka  +  ck3  +  dk  (2-13) 

where  b.  ,  c,  ,  and  d.  are  real  constants. 
k   k       k 

Substitution  of  equation  (2-13)  into  equation  (2-12)  yields: 

(2-14) 


Sua  +  c.3  +  D,  =  0 


where 


B„a  +  C2B  +  D„  =  0 


Bl  ■  j0  (-1)kVnUk-l(5) 


B2  =  X  (-1)kVnUk(«» 

k=0 


k=0 


n      ,    , 
C2  =  I       (-l)Kck03^Uk(c)  (2-15) 

k=0 


11 


Dl  =  Jn  <-1)kVnUk-l(^ 
k=0 

D2  =  J.  '-"Vn^1^    •  (2-15> 

k=0 

Solution  of  equation  (2-14)  by  Cramer's  Rule  yields  the 

following  parameter  plane  solution  equations: 

C1D2-C2D1 

a  = 


B1C2"B2C1 


B  m   B2D1-B1D2 
B1C2"B2C1 


(2-16) 


The 


functional  dependence  of  B, ,  B~ f    C, ,  C~,  D.. ,  D~  on 


£  and  co  will  be  omitted  for  simplicity  in  this  report.   A 
list  of  the  algebraic  form  of  the  Chebyshev  functions  of  the 
first  kind  U,(c;)  is  given  in  Appendix  I.   Appendix  II  gives 
a  table  of  U,  (C)  values  for  selected  values  of  ?. 

The  parameter  plane,  as  defined  by  Siljak  (2)  is  a  rec- 
tangular coordinate  plot  with  a  as  the  abscissa  and  3  as  the 
ordinate.   Equation  (2-14)  permits  the  mapping  of  points, 
excepting  real  axis  points,  from  the  s-plane  into  the  param- 
eter plane. 

For  example,  £  may  be  fixed  £  =  £  ,  and  varying  co   from 
zero  to  infinity  results  in  a  constant  zeta  curve  in  the 
parameter  plane.   This  curve  specifies  the  a, 3  pairs  which 
will  cause  equation  (2-1)  to  have  a  pair  of  complex  roots 
with  damping  ratio  £  =  C _•   Similarly  fixing  co   =  co    and 

O  XI         X 1  S3 

varying  c  from  zero  to  one  will  produce  a  curve  on  the  param- 
eter plane  which  will  specify  the  a, 3  pairs  for  which  a  pair 
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of  complex  roots  in  equation  (2-1)  have  a  constant  natural 

frequency,  co  .   It  should  be  noted  that  this  mapping  is  not 

conformal  since  a  pair  of  complex  roots  on  the  s-plane  maps 

into  one  point  on  the  parameter  plane  and,  as  will  be  seen, 

a  real  axis  point  on  the  s-plane  maps  into  a  straight  line 

on  the  parameter  plane. 

For  real  axis  points  in  the  s-plane  note  that  ?  =  1  and 

letting  oo   =  -a ,  the  equation  for  s  becomes  s  =  -a.   The 

characteristic  equation  becomes: 

n        , 
F(s)  =  I      a,  (-a)K  =  0   .  (2-17) 

k=0   K 

Substitution  of  equation  (2-13)  into  equation  (2-17)  yields: 

n  , 

•  F(s)  =  I       (b.  a+c,  3+d,  )  (-a)K  =  0   .       (2-18) 
k=0   K   K   K 

Simplifying  equation  (2-18)  yields: 

aB(a)  +  $C(o)    +   D(a)  =  0  (2-19) 


where : 


n        k    k 

B(a)  =  I       (-l)Kb,a 
k=0       K 

5     k   k 

C(a)  =  I      (-l)Kc,aK  (2-20) 

kio  k 

D(a)  =  I       (-l)kd,ak 
k=0       K 

For  a  fixed  value  of  a,  B(a)  ,  C(a)  ,  and  D(a)  are  constant  and 
equation  (2-19)  in  the  equation  of  a  straight  line  on  the  pa- 
rameter plane.   The  a, 3  pairs  corresponding  to  this  line  cause 
equation  (2-1)  to  have  a  real  root  at  s  =  -a. 
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B.   MAPPING  OF  THE  PARAMETER  PLANE 

The  mapping  of  the  s-plane  into  the  parameter  plane  per- 
mits a  designer  to  choose  or  adjust  parameters  so  that  charac- 
teristic equation  roots  lie  within  specified  areas  of  the 
s-plane. 

The  rules  and  graphical  techniques  for  the  mapping  proc- 
ess are  discussed  in  detail  by  Siljak  (2)  and  Thaler  (13) . 

The  major  mapping  rules  and  interpretation  of  the  result- 
ing curves  are  illustrated  in  the  following  example. 

Consider  a  linear  fourth  order  system  with  the  follow- 
ing characteristic  equation: 

F(s)  =  s4  +  (ct+15)s3  +  (15a  +  50)s2  +  (50a)s  +  6006  =  0  . 
Figure  II-l  shows  the  mapping  of  the  characteristic  equa- 
tion in  the  parameter  plane.   Inspection  of  Figure  II-l  re- 
sults in  the  following  observations : 

1)  A  pair  of  complex  conjugate  points  on  the  s- 
plane  maps  into  a  single  point  on  the  param- 
eter plane. 

2)  A  real  axis  point  on  the  s-plane  maps  into  a 
straight  line  on  the  parameter  plane. 

3)  The  stable  area  on  the  parameter  plane,  for 
this  polynomial  is  defined  by  the  positive  a 
axis,  the  zeta  equal  zero  curve,  and  infinity. 

4)  Complex  roots  of  the  polynomial,  for  specific 
values  of  a  and  3  are  located  in  the  region 
on  the  parameter  plane  defined  by  the  zeta 
equal  zero  and  zeta  equal  one  curves. 
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For  example,  point  M.  is  at  the  intersection 

of  the  co   =2.0  and  C  =  .4  curves  therefore 
n 

the  characteristic  equation  has  roots  at: 


s  =  -(.4)  (2.0)  ±  j2.0/l-(.4)  2 

=  -0. 8  ±  jl.83   . 
5)   The  remaining  portion  of  the  stable  area 
defined  by  the  zeta  equal  one  curve,  the 
positive  a  axis,  and  infinity  is  a  real 
root  area.   For  example  point  M,  is  at  the 
intersection  of  four  real  root  lines  giving 
roots  of: 

s  =  -1.0,  -2.0,  -8.0,  and  -9.0 
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x  scale  =  l.o   units/in 
V  scale  =  .10   units/in 

s  (a+15)s3  +  (15a+50)s2  +  50as  +  600S  =  0 
Figure  Il-l    Parameter  Plane  Example 
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III.   SINGULAR  LINES  ON  THE  PARAMETER  PLANE 

The  parameter  plane  developed  by  Siljak  (2)  and  the  com- 
puter program  written  by  Nutting  (7)  solved  for  all  of  the 
roots  of  the  characteristic  equation  according  to  then  exist- 
ing theory.   However,  Bowie  (8)  showed  that  for  some  systems 
the  existing  parameter  plane  method  did  not  solve  for  all  of 
the  roots  of  the  characteristic  equation.   Bowie  found  that 
these  undetermined  roots  were  actually  located  on  a  constant 
zeta-constant  natural  frequency  line  which  has  been  defined 
as  a  singular  line. 

A.   MATHEMATICAL  BASIS  OF  SINGULAR  LINES 

As  pointed  out  in  Section  II  the  solution  of  the  parameter 
plane  values  (equation  2-16)  was  carried  out  by  the  applica- 
tion of  Cramer's  Rule  to  equation  (2-14).   The  application  of 
Cramer's  Rule  to  equation  (2-14)  is  valid  only  when  the  aug- 
mented coefficient  matrix  of  this  equation  is  non  singular  or 
of  rank  equal  two.   Bowie  showed  that  when  the  augmented 
coefficient  matrix  of  equation  (2-14)  is  singular  and  the 
rank  is  equal  to  one,  a  condition  exists  such  that  8  becomes 
a  linear  function  of  a  and  a  singular  line  results.   In  or- 
der for  the  augmented  coefficient  matrix  to  be  singular  and 
of  rank  equal  one  the  following  conditions  must  prevail: 
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B. 


B, 


-D. 


-D. 


B1C2  "  B2C1  =  ° 


C1D2  "  C2D1  =  ° 


(3-1) 


(3-2) 


Bl   "Dl 


B2   -D2 


■  B2D1  "  B1D2  =  ° 


(3-3) 


For  example,  let  the  solution  of  equations  (2-15)  qive 
the  following  values  for  a  specific  r, ,    to   pair: 


Bi  =  2 


C.  =  5 


Dl  =  6 


B2  =   4 


c2  -  10 


D2  =  12 


Equations  (3-1,  3-3)  are  satisfied  and  equations  (2-14) 
become : 

2a  +   53  +   6=0 

4a  +  103  +  12  =  0   .  (3-4) 

The  resulting  relationship  between  a  and  3  is: 

Q      2a+6 

p  =  -  — ? 


which  is  a  straight  line  on  the  parameter  plane.   This  line 
is  a  locus  of  a, 3  pairs  which  will  cause  the  system  charac- 
teristic equation  to  have  a  pair  of  fixed  complex  roots  with 
values  of  £  and  oo   corresponding  to  those  used  in  solvinq 
equation  (2-15)  above. 


B.   METHODS  FOR  SHOWING  EXISTENCE  OF  SINGULAR  LINES 

In  general  one  or  more  singular  lines  may  be  shown  to  ex- 
ist for  a  particular  system  by  the  application  of  equations 
(3-1,  3-3)  however,  this  approach  requires  the  substitution 
of  specific  £,  co   pairs  into  equation  (2-15).   Since  there 
are  an  infinite  number  of  C ,  co   pairs  which  might  satisfy 
equations  (3-1,  3-3),  it  is  apparent  that  additional  criteria 
for  the  possible  C ,  co   pairs  must  be  established. 

Method  I 

Bowie  derived  the  following  method  for  finding  singular 
lines  on  the  parameter  plane  (8).   Proof  of  the  following 
equations  will  not  be  shown  but  may  be  found  in  reference  (8). 

The  expansion  of  B.C2  -  B-C,  is  shown  as  equation  (3-5) 
on  the  following  page,  or  in  compact  notation 

n 

B,C0  -  B0C,  =   T   ce.   =  0  (3-6) 

12     2  1    i^1    ij 

j  =  l 

where : 

CE  =  [ce. .]  is  an  (nxn) 

upper  trangular  array  whose  elements  are: 

[b.  ,c.  .  ,-b.  .  nc.  J  (-1)1"1U.  U)a)2l+i~2   for  all  i+j  <  n+2 
l-l  i+j-1   i+j-1  i-I         j  n  J 


for  all  i+j  >_  n+2 
(3-7) 


Similarly : 

n 


C,D0  -  C0Dn  =  I      ce. .  =  0  (3-8) 

12     2  1    _.  *- .,    i  j 


i=l 

j  =  l 
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where : 


CE  =  [ce. .]  is  an  (nxn) 


upper  triangular  array  whose  elements  are: 

[c.  ,d.  .  ,-c.,.  ,d.  .] (-1)1"1U. (c)w2l+j"2   for  all  i+j  <  n+2 
l-l  i+n-1   l+i-l  l-l         i     n  J 


and 


for  all  i+j  >_  n+2 
(3-9) 


n 
B1D2  ~  B2D"'  =  "  ^   ce- -  =  °  (3-10) 


!D1  "  "2   CGij 
1=1 

j  =  l 


where : 


CE  =  [ce. .]  is  an  (nxn) 
upper  triangular  array  whose  elements  are 

[b.  nd.  .  ,-b.L.  .d.  .  ]  (-1)1_1U.  (c)u)  2l+J~2   for  all  i+j  <n+2 
l-l  l+i-l   l+i-l  i-I         3  n  J 

0  for  all  i+j  _>  n+2 

(3-11) 

Equations  (3-6),  (3-8)  and  (3-10)  give  three  polynomials 

in  go   and  C  of  order  2N-1  in  co   and  N-l  in  C.   For  simple  sys- 

n      /  n 

terns  it  is  possible  to  solve  for  unique  values  of  £  and  to 
which  satisfy  equations  (3-6) ,  (3-8)  and  (3-10) .   For  more 
complex  systems  it  is  necessary  to  set  £  =  £   in  one  of  the 
equations ,  solve  for  to  ,  and  check  the  resulting  values  in 
the  two  remaining  equations.   This  approach  is  best  handled 
in  a  digital  computer. 
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Example  I 

Let  the  system  characteristic  equation  be: 

F(s)  =  s3  +  as   +  3s  +  2.5g  -  12.5  =  0  . 
Solving  equations  (3-6),  (3-8),  and  (3-10)  the  follow- 
ing is  obtained: 

B1C2  "  B2C1  =  5?%  "  wn  =  ° 

5  2 

2    1  12  n  n 

C,D0    -    C0D.    =    -2w4C    +    2.5o)3(4c2-D    +    12  .  5oo      =    0       . 
12  21  n^  n  n 

Solving  the  above  for  £  and  co   the  following  unique  val- 


ues are  obtained: 


5  =  —   co   =  /5 


/5     n 

Applying  these  values  to  equation  (2-14)  results  in: 
5a  -  2.56  =  -2.5 
and 

2/5  a  -  /5  3  =  -/5 
and  the  singular  line  relation  is: 

3  =  2a  +  1 

Method  II 

The  following  approach  to  the  singular  line  problem  was 
developed  by  Cadena  (10) . 

Assuming  that  a  system  has  a  singular  line  for  some  £, 
co   pair,  there  must  be  a  pair  of  complex  roots  which  are 
fixed  for  all  a, 3  pairs  in  the  parameter  plane.   Therefore, 
let  a  pair  of  complex  roots  of  the  characteristic  equation 
be  given  by : 
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s   =    -a    ±    j/oo2-a2  (3-12) 

J      n 

where   a    =    ?a3    .      After  multiplication   and  manipulation   the 

equation    (3.12)    becomes: 

s2    +    2a   +    co2   =    0    .  (3-13) 

n 

If  this  pair  of  roots  is  divided  into  the  characteristic 
equation  there  results  a  quotient  which  is  in  effect  a  re- 
duced characteristic  equation  which  has  as  roots  the  remain- 
ing roots  of  the  original  characteristic  equation.   Also 
there  is,  in  general,  a  remainder  in  the  form: 

I      sK[f  (a,3,a,co  )]   .  (3-14) 

k=0  n 

If  the  complex  roots,  chosen  as  the  divisor,  are  to  be 
roots  of  the  original  characteristic  equation,  the  coeffi- 
cients of  equation  (3-14)  must  be  independently  equal  to 
zero.   Setting  these  coefficients  equal  to  zero  with  the  re- 
striction that  |3  is  a  linear  function  of  a  results  in  a 
system  of  polynomials  which  is  the  same  as  the  system  devel- 
oped in  Method  I  above. 

Example  II 

Dividing  the  characteristic  equation: 

F(s)  =  s3  +  as2  +  3s  +  2.53  -  12.5  =  0 

(from  Example  I  above)  by: 

2  2 

s   +  2as  +  03 
n 

results  in  a  quotient  or  reduced  characteristic  equation  of: 

s  +  a  -  2  =  0 


23 


and   a   remainder  of 

(3-oo2-2aa+4a2)  s   +    (2 .  53-aoo2+2aoo   -12.5)    =    0    . 
n  n  n 

Since  both  terms  of  the  remainder  must  be  equal  to  zero 

and  since  3  is  a  linear  function  of  a  the  derivative  of  3  with 

respect  to  a  must  be  the  same  for  both  terms.   From  the  first 

term: 

da 
and  from  the  second  term 

/ID       2(jJ 

d3     n 


da     5 

Equating  the  derivatives : 

2 

oo   =  5a 
n 

Substitution   of   a    =    Coo      yields: 

n  •* 

oj   =  5£ 
n 

Returning  to  the  first  remainder  term,  solving  for  3  and  sub- 
stituting for  3  in  the  second  remainder  term  results  in: 

a(10a-2oo2)  +  5co2  -  20a2  -  25  +  4aoo2  =  0   . 
n      n  n 

Since  the  above  equation  must  be  satisfied  for  all  a  it  must 

be  true  that  the  a  coefficient  is  equal  to  zero  and  that  the 

remaining  constant  terms  sum  to  zero.   From  the  a  coefficient: 

2 

oj  =  5a 
n 

which  was    derived   above    using   derivatives. 

Summing   the    constant    terms : 

5oo2    -    20a2    -    25   +    4aoo2    =    0    . 
n  n 
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Returning  to  the  first  remainder  term,  solving  for  a  in 
terms  of  3,  substituting  into  the  second  remainder  term  and 
setting  the  sum  of  the  constant  terms  equal  to  zero  results  in 

4     or 

co   =  25a 
n 

Summarizing  the  above  relationships  and  letting  a  =  £co 

results  in: 

to   =  5C 
n 

un  =  I?  [5co^(4c:2-l)+25] 

3 

co   =  25c 
n 

Comparing  these  results  with  the  results  of  Example  I 

above  shows  that  the  singular  line  relationships  are  identical 

Solving  the  polynomials  for  Z,    and  co   results  in: 

C-V 

co   =  /5 
n 

a   =    l 

Replacing  these  values  in  one  remainder  term  provides 

the  required  a,  3  relationship  for  the  singular  line: 

3  =  2a  +  1   . 

The  third  root  is  obtained  by  substitution  of  a  into  the 

quotient  term  giving: 

s  =  -a  +  2 

Method  III 
Karmarker  derived  the  following  method  for  showing  the 
existence  of  singular  lines  in  the  parameter  plane  (11) . 
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The  system  characteristic  equation  may  be  written  as: 

(3-15) 


n       n-1 
s   +  a   s 
n 


+  a,  s  +  a  =0 
1     o 


where : 


ak  "  bka  +  ck2  +  dk 


(3-16) 
Let  one  pair  of  complex  roots  for  equation  (3-15)  be: 


s  =  -a  ±  j/co  -a: 
J   n 


(3-17) 
Now  equation  (3-15)  may  be  factored  into  two  polynomials: 

(3-18) 
which  is  obtained  from  equation  (3-17)  and: 

=  0 


2  2 

s   +  2as  +  o)   =  0 
n 


n-2       n-3 
s    +  x,  s 


x   _s  +  x   0 

n-3     n-2 


(3-19 


which  is  the  reduced  characteristic  equation. 

Multiplying  equation  (3-18  and  (3-19)  equating  to  the 
coefficient  in  equation  (3-15)  and  putting  the  matrix  nota- 
tion results  in: 

[A]  [B]  =  [C] 
where 


A  = 


b    1 

n-1 

c 

n- 

-1 

-1 

0 

0       .     . 

0 

bn      9 

n-z 

c 
n- 

-2 

-2a 

-1 

0 

• 

• 

• 

2 

-co 
n 

-2a 

-1 

• 

• 

• 

0 

2 

-co 
n 

-2a 

• 

• 

• 

• 

0 

2 

-co 
n 

• 

• 

« 

• 

• 

0 

-1 

bl 

cl 

0 

0 

0 

-2a 

b„ 

c« 

0 

0 

0 

2 

-co 

n 


(3-21) 
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B  = 


a 
6 

x] 


X 


n-3 


x 


n-2 


C  = 


2a  -  d 


n-1 


w   -  d 
n    n-2 


-d 


n-3 


-d. 


-d, 


(3-21) 


For  singular  lines  to  exist  3  must  be  a  linear  function 
of  a.   This  criterion  plus  the  techniques  of  linear  algebra 
result  in  the  following  singular  line  existence  theorem: 

Singular  Line  Existence  Theorem 

A.  All  determinants  of  order  n  that  can  be  formed 
from  the  augmented  matrix  [A,C]  are  zero;  (n+1)  such 
determinants  may  be  formed.   Note  that  n  of  these 
form  an  independent  set. 

B.  At  least  one  of  the  determinants  of  order  (n-1)  , 
that  can  be  formed  from  the  augmented  matrix  [A,C] 
must  be  non-zero.   There  are 

n+1      n 
n-1     n-1 
such  determinants. 
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The  conditions  that  must  exist  for  singular  lines  may  be 
found  using  the  existence  theorem  as  shown  in  the  following 
example. 

Example  III 

Consider  the  characteristic  equation  used  in  Example  II: 

F(s)  =  s3  +  as   +  3s  +  2.56  -  12.5  =  0 


[A]  = 


1 

0 

-1 

0 

1 

-2a 

0 

2.5 

2 

-oo 

n 


[B]  = 


[C]  = 


2a 
2 

00 

n 
12.5 


The  (n+1)  determinants  of  order  n  which  can  be  formed 
from  the  augmented  matrix  [A,C]  are: 


1  0 
0  1 
0   2.5 


-1 

-2a 
2 


-oo 


n 


2a 
2 


oo 


n 


0   2.5   12.5 


1 

2a 

-1 

0 

2 

00 

n 

-2a 

0 

2.5 

2 

-00 

n 


2a 
2 


CO 


n 


-1 

-2a 
2 


12.  5  2.  5   -oo 


n 


Solving  the  above  determinants  and  equating  each  to  zero 
results  in  the  following: 


oo  =  5a 
n 


oo   =  25a 
n 
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2         10a2+12.5 
n   '        2a+2 . 5 


2         _ 
co      =    5 
n 


Solving  for  a  and  co  : 

3  n 

a  =  1 

co   =  /5 
n 

It  should  be  noted  that  replacing  a  with  Ceo   in  the  first 
three  equations  given  above  results  in  the  same  set  of  equa- 
tions as  given  in  Examples  I  and  II.   The  last  equation  is  re- 
dundant. 

Substituting  for  a  and  co   in  equation  (3-20)  and  solving 
results  in: 

6  =  2a  +  1 
r~  =  x  =  a  -  2 

Again,  this  is  the  same  result  obtained  in  Examples  I 
and  II  above. 

Method  IV 
The  following  parameter  plane  equation  and  singular  line 
theory  was  derived  by  Randolph  (12). 

A  second  derivation  of  the  parameter  plane  transformation 
which  is  somewhat  different  than  that  described  in  Section  II 
is  described  below. 

Let  the  system  characteristic  equation  be  written  as : 

n 

r 

L 
k=0 


F(s)  =  t     aksk  =  °  (3-22) 


29 


where 


ak  =  bka  +  cke  +  dk   '  (3-23) 

The  complex  variable  s   may  be  written  as : 

s   =  co  e-1    =  co  [cos  ke+jsin  k0]  (3-24) 

n        n 

where  6  is  the  angle  from  the  positive  real  axis  to  a  Doint 

in  the  s-plane  and  co   is  the  natural  frequency  of  the  system. 

Substituting  equation  (3-24  and  (3-23)  into  equation 

(3-22)  results  in: 

5  k 

F(s)  =  I       [b,  a+aS+d,  ]co  [cos  ke+jsin  k9]  =  0   .  (3-25) 
k=0    K  n 

Letting  the  real  and  imaginary  terms  go  independently  to 

zero  results  in 

Bxa  +  C^B  +  Dx  =  0 

B2a  +  C2e  +  D2  =  0  (3-26) 

where 

n     ,  n     . 

B,  =   T   b,  co   cos  k9  B9  =   7   b,  co   sin  k9 
1    k=0  K   n  *        k=0   K  n 

r     k  ?  k 

C,  =  I      c,co   cos  k0         C9  =  I      c,co   sin  k0 

x   k=0     n  z    k=0   K  n 

n     .  n     , 

D   =   7   d,  co   cos  k9         D0  =   7   d,  co   sin  k6  .  (3-27) 
1    k=0   K  n  Z        k=0   K  n 

It  should  be  noted  that  the  coefficients  of  equation 
(3-26)  are  not  the  same  as  described  in  Section  II  above. 

From  the  theory  of  linear  algebra  it  may  be  shown  that 
two  functions,  f-jta,^)  and  f2(a,3)  are  linearly  dependent  if 
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and  only  if  K^^S)  +  K1f1(a,3)  +  K2f2(a,3)  =  0  for  all 
a,  3  in  a  domain  common  to  f,  and  f~  for  some  K..  and  K~  not 
both  zero. 
Let: 

f1  (a, 6)  =  B.,a  +  C,3  +  D1 

f2(a,8)  =  B2a  +  C23  +  D2   .  (3-28) 

Then  a  necessary  and  sufficient  condition  that  f, (a, 3) 
and  f?(a/3)  be  linearly  dependent  is  that  there  exists  a  con- 
stant M  such  that: 


B2 

=    MB. 

C2 

=   MC, 

D2 

=   MC, 

(3-29) 

The  conditions  that  must  be  satisfied  in  order  for  a 
singular  line  to  exist,  may  be  obtained  as  follows: 

(a)  Write  equation  (3-27)  . 

(b)  Relate  these  equations  using  equations  (3-29)  . 

(c)  Solve  for  M,  9,  and  go 

(d)  The  resulting  singular  line  equation  is  obtained 
from  equation  (3-2  6) 

Example  IV 

Using  the  same  characteristic  equation  as  in  previous 
examples : 

F(s)  =  s3  +  as2  +  3s  +  2.53  -  12.5  =  0 
Equations  (3-29)  may  be  written  as 

MB,    =   M[to      cos29]    =    oj2    sin20    =    B„ 
In  n  2. 
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MC,    =   M[2.5+co      cosG]    =    co      sine    =   C~ 
l  n  n  z 

MDn    =   M[-12.5+co3    cos36]    =    oj3    sin39    =   D„ 
1  n  n  2 

Solving    for  M,    9,    and   co      results    in: 

M   =    i 
3 

6   =    cos_1(-   — ) 


/5 


co      =    /5 
n 


Since : 


=   cos       (-C) 


This    results    in   a   complex  pair  of   roots    at: 

s   =    -1    ±    j2 

Returning  to  equation  (3-26)  and  solving  for  the  singu- 
lar line  equation  results  in: 

3  =  a  -  2 

Dividing  the  characteristic  equation  by  the  complex 
roots  results  in  the  third  root  of  the  system: 

r  =  a  -  2 

These  results  are  the  same  as  those  of  previous  examples. 

The  preceding  portion  of  this  section  developed  methods 
for  showing  the  existence  of  singular  lines  and  the  illus- 
trative examples  indicated  procedures  for  finding  the  neces- 
sary conditions  for  singular  lines  using  the  existence  the- 
orems.  These  procedures  work  well  for  simple  systems,  but 
become  complicated  for  the  higher  order  systems  since,  in 
general,  the  order  of  con  is  2  (n-1)  and  the  order  of  a  is 
n-1  in  the  resulting  equations. 
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The  following  procedure  provides  a  simple  technique  for 
finding  singular  line  characteristic  equations.   Manipulation 
of  the  resulting  characteristic  equations  will  in  general  pro- 
vide a  singular  line  plant. 

Let  a  system  characteristic  equation  be  written  as: 

n 

1 
k=0 

where : 


F(s)  =  I      skak  =  0  (3-30) 


ak  =  abk  +  eck  +  dk   '  (3-31) 

Assume  that  F(s)  is  a  singular  line  characteristic 

equation  for  some  £  =  £  ,  co   =  go   .   Then  equation  (3-30) 

may  be  rewritten  in  factored  form  as : 

F(s)  =  (s2+2c  co   +co2  )  [  (s+r,  )  (s+r0) (s+r    )]  =  0 

s  ns   ns       1      Z  n-2 

(3-32) 
In  general  the  roots  r, ,  r~ ,  ....  r  _2,  are  functions 
of  a  and  3.   However  for  the  singular  line  case  a  and  3 
must  appear  in  linear  terms  such  that  the  order  of  a  and  3 
is  one  and  a3  products  do  not  appear.   Therefore  only  one 
root  of  r,r2....r  _2,    may  contain  terms  in  a  and  3  and  this 
root  must  be  real.   The  remaining  roots  may  be  fixed  at  any 
desired  value.   If  the  resulting  mobile  real  root  contains 
only  a  as  a  variable  a,  will  be  a  function  of  a.   A  singular 
line  may  now  be  obtained  by  substituting  3  as  a  linear  func- 
tion of  a  into  one  or  more  of  the  resulting  a,  terms. 

For  example  suppose  that  for  a  third  order  system  a 
singular  line  is  desired  such  that  the  roots  are: 
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s  =  -1  ±  j2 

s  =   a  -  2 
Substituting  into  equation  (3-32)  yields: 
F(s)  =  (s  +2s+5)  (s+a-2)  =  s3  +  as2  +  (2a+l)s  +  5a  -  10  =  0  . 

(3-33) 
Let  3  =  2a  +  1  and  substitution  into  equation  (3-33)  results 
in: 

F(s)  =  s3  +  as2  +  3s  +  2.53  -  12.5  =  0   .     (3-34) 
Equation  (3-34)  is  a  singular  line  characteristic  equa- 
tion for  3  =  2a  +  1  as  shown  in  previous  examples. 

The  procedure  indicated  above  results  in  a  unique  sin- 
gular line  for  a  given  characteristic  equation.   A  more  gen- 
eral case  is  examined  by  Bowie  (8)  in  which  there  are, 
apparently,  an  infinite  number  of  singular  lines  for  a  given 
characteristic  equation.   The  system  examined  is  a  complex 
cross-coupled  type  with  powers  of  s  to  the  sixth  power  occur- 
ring in  the  characteristic  equation.   There  does  not  appear 
to  be  a  tractable  method  for  producing  similar  systems  at 
the  present  time. 
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IV.   EXAMPLES  OF  SYSTEMS  WITH  SINGULAR  LINES 

This  section  contains  block  diaqrams,  characteristic 
curves,  and  step  response  characteristics  for  systems  which 
have  singular  lines.   The  systems  are  grouped  by  the  nature 
of  their  characteristic  equations. 

A.   SYSTEMS  WITH  ONE  SINGULAR  LINE  AND  ONE  MOBILE  ROOT 


Characteristic 

Sinqular 

Ex. 

Equation 

Roots 

Line 

Figs . 

1 

s3+as  +(2a+l)s+6=0 

s=-l±j2 
s=- (a-2) 

8=5(a-2) 

IV-1 
IV- 2 

2 

s  +as  +  6s+2 . 56-10  0= 

=  0 

s=-l±j2 
s=-(a-2) 

6=2a+l 

IV- 3 
IV- 4 

s3+as2+(43+5)s+3=0      s=-l±j2  B=5(a-2)     IV-5 

s=-(a-2)  IV-6 

s4+(a+10)s3+(12a+l)s2    s=-l±j2  6=5(a-2)     IV-7 
+  25as+10B=0     s=-(a-2)  IV- 8 

s=-10 

s4+(a+10)s3+(12a+l)s2    s=-l±j2  6=25a       IV-9 
+  6s+26-100=0    s=-(a-2)  IV-10 

s=-10 

s4+(a+10)s3+(12a+l)s2    s=-l±j2  6=5(a-2)     IV-11 
+  (56+50) s+10B=0     s=-(a-2)  IV-12 

s=-10 


As  a  is  increased  the  mobile  root  moves  further  from  the 
origin  of  the  s-plane  and  the  effect  of  this  root  becomes 
negligible  for  values  of  a  greater  than  20.   This  trend  is 
shown  by  the  step  response  characteristics  for  the  respective 
systems . 
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The  movement  of  the  mobile  root  along  the  real  axis  of 
the  s-plane  maps  along  the  singular  line  in  the  parameter 
plane;  this  is  shown  by  substituting  x  =  -(a-2)  into  the  sys- 
tem characteristic  equation.   In  the  fourth  order  systems, 
the  fixed  real  root,  s  =  -10,  also  maps  along  the  singular 
line. 
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x   scale 
y    scale 


002  004 
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s3+as2+(2a+l)s+6=0 


Figure  IV-2   Parameter  Plane  Curves  -  Example  1 
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x  scale  =  20  units/in 
y  scale  =  20  units/in 

s3+as2+Bs+2. 58-100=0 
Figure  IV-4   Parameter  Plane  Curves  -  Example  2 
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x  scale 
y  scale 


20  units/in 
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s3+as2+(46+5)s+6=0 


Figure  IV-6   Parameter  Plane  Curves  -  Example  3 
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Figure  IV-8   Parameter  Plane  Curves  -  Example  4 
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Figure  IV-14   Parameter  Plane  Curves  -  Example  7 
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Figure  IV-16   Parameter  Plane  Curves  -  Example  8 
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C.   SYSTEM  WITH  LOCKED  ROOTS 

The  system  shown  in  Figure  IV-17  has  the  following  char- 
acteristic equation: 

s3  +  (a+453+22)s2  +  (22a+990B+45) s  +  100  =  0  . 
Setting  a  =  -453  forces  the  system  to  have  the  following 
roots  : 

s  =  -20 
s  =  -1  ±  j2 
There  are  no  parameter  plane  curves  for  this  system  since 
the  parameter  consists  of  an  infinite  number  of  singular 
lines . 
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V.   SELF-ADAPTIVE  CONTROL  USING  SINGULAR  LINES 

One  important  application  of  singular  line  systems  is 
that  of  self-adaptive  control.   For  example,  a  parameter 
which  is  a  function  of  environmental  conditions  could  be 
compensated  for  by  linearly  varying  a  second  parameter. 
This  situation  might  occur  in  an  aircraft  where  the  dynamic 
flight  characteristics  are  a  function  of  altitude. 

Consider  Figure  V-l,  below,  which  shows  a  singular  line 

for  L   =    C   and  00=00   on  the  parameter  plane.   3  is  the 
s      n    ns        "         * 

varying  plane  parameter  and  a  is  an  adjustable  parameter  in 
the  controller.   The  initial  operating  condition  is  on  the 
singular  line  at  point  P,  so  that  dominant  complex  roots 
are  at 


s  =  -c  u   ±  joo   /l-c;2 
s  ns    J  ns    s 

When  3  is  varied  by  an  amount  A3  the  operating  point  is 

changed  to  point  P2  by  changing  a  by  an  amount,  Aa.   The  sys' 

tern  maintains  roots  at 


s  =  -C  00    ±  joo   /1-C2 
s  ns    J  ns     s 

The  other  characteristic  roots  may  be  changed,  but  for  rea- 
sonable changes  in  a   and  3  root  dominance  is  retained. 


3 


singular  line 

£  =  £    00   =00 

s    n    ns 


^  a 


Figure  V-l   Basic  Concept  of  Self-Adaptive  Control 
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The  system  in  Figure  V-2  is  an  example  of  a  singular  line 
plant  (pole-zero  cancellation)  which  has  been  designed  for 
self-adaptive  control  by  a  procedure  developed  by  R.  A. 
Desrosiers  (13) .   Explanation  of  the  various  control  func- 
tions will  not  be  explained  but  may  be  found  in  reference 
(13) .   Figure  V-2  also  shows  the  step  response  of  the  system 
with  a  12%  step  change  of  3  at  t  -  .09  sec.   As  indicated, 
the  adapted  plant  and  model  responses  coincide.   (The  model 
is  the  basic  plant  with  no  parameter  changes.) 
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Figure  V-2   Step  Responses  for  Self-Adapting  Plant 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  WORK 

Feedback  control  systems  containing  singular  lines,  while 
occurring  in  a  limited  number  of  cases ,  provide  a  tool  for 
the  design  engineer  working  in  self-adaptive  control  and 
other  areas.   Methods  for  predicting  the  existence  of  singu- 
lar lines  have  been  shown  with  examples  given  in  the  use  of 
these  methods.   As  indicated  in  previous  sections  of  this 
report,  the  equations  resulting  from  the  application  of  known 
existence  techniques  are  complex  and  do  not,  in  general,  pro- 
vide a  tractable  means  of  finding  singular  line  systems. 

From  the  standpoint  of  the  design  engineer,  it  would  be 
desirable  to  formulate  a  general  set  of  rules  which  would 
provide  a  specific  system  with  singular  lines.   Research  to 
date  does  not  indicate  the  existence  of  such  a  general  set 
of  rules.   Further  development  is  needed  in  this  area. 

A  second  approach  might  be  to  design  a  number  of  singu- 
lar line  systems  using  characteristic  equations  as  starting 
points.   A  set  of  typical  singular  line  systems  would  then 
be  made  available  to  the  design  engineer  for  correlation  be- 
tween a  specific  system  and  possible  singular  line  charac- 
teristics . 

Another  area  for  future  work  might  be  to  develop  further 
theorems  for  singular  lines  on  the  parameter  plane.   For  ex- 
ample, all  real  root  lines  for  singular  line  systems  were 
found  to  lie  on  the  singular  line.   Further  study  may  Drove 
that  this  is  true  in  general. 
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APPENDIX  I 

TABLE  OF  THE  CHEBYSHEV  FUNCTIONS  OF  THE  SECOND  KIND 

The  Chebyshev  functions  of  the  second  kind,  U,  (£),  are  de- 
fined by  the  recursion  relation: 


where 


Uk+lU)  "  2cUkU)  +  Uk-lU)  =  ° 


uQ(c)  =  o 


u1(c)  =  1  . 

These  functions  for  k  =  -1,...,  10  are 

k  Uk(^ 


-1  -1 

0  0 

1  1 

2  2c 

3  4?2  -  1 

4  8C3  -  4c 

5  16?4  -  12c2  +  1 

6  32£5  -  32C3  +  6c 

7  64C6  -  80C4  +  24C2  -  1 

8  128C7  -  192C5  +  80C3  -  8c 

9  256C8  -  448C6  +  240c4  -  40C2  +  1 

10  512c9  -  1024c7  +  672c5  -  160c3  +  10c 
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